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We report a search for the anomalous production of events with multiple charged leptons in pp col- 
lisions at t/s — 1.96 TeV using a data sample corresponding to an integrated luminosity of 346 pb _1 
collected by the CDF II detector at the Fermilab Tevatron. The search is divided into three-lepton 
and four-or-more-lepton data samples. We observe six events in the three-lepton sample and zero 
events in the >4-lepton sample. Both numbers of events are consistent with standard model back- 
ground expectations. Within the framework of an R-parity violating supergravity model, the results 
are interpreted as mass limits on the lightest neutralino (x?) and chargino (xf) particles. For one 
particular choice of model parameters, the limits are M(xi)>110 GeV/c 2 and M(xf)>203 GeV/c 2 
at 95% confidence level; the variation of these mass limits with model parameters is presented. 
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Numerous attempts to resolve theoretical problems 
with the standard model of elementary particle physics 
(SM) require the existence of new particles with masses 
at the electroweak scale, ~100 GeV/c 2 HHHQ. At the 
Tcvatron and LHC hadron colliders, the production cross 
sections for these particles are predicted to be orders 
of magnitude smaller than for SM processes. Analysis 
methods for reducing background levels while preserving 
new physics signals include searching for leptons (£) [5| 
with large momentum transverse to the beam axis (pr), 
like-sign leptons, or signatures with transverse energy im- 
balance ($t)- Requiring several leptons effectively re- 
duces the SM backgrounds while maintaining sensitivity 
to low i^x regions. In this Letter, we present a search 
for an excess of events containing three or more leptons 
above the SM prediction. We use a data sample with 
an integrated luminosity of 346 pb" 1 of pp collisions at 
^/s = 1.96 TeV collected by the CDF II detector from 
March 2002 to August 2004. 

Several new physics models predict final states with 
four or more leptons that can be produced in pp collisions, 
including R-parity violating supersymmetry (i£ p SUSY) 
non- minimal supersymmetric models (nMSSM) Q, 
and production of doubly-charged Higgs pairs Q. In 
ftp SUSY, the lightest supersymmetric particle (LSP) 
may decay directly into two leptons and a neutrino. With 
a pair of LSPs in each event, there could be at least 
four leptons plus j£r- Another possibility is the nMSSM, 
which includes an expanded Higgs sector and their super- 
symmetric partners. "Cascade" decays of these particles 
may result in multiple leptons, jets, and less than in 
other SUSY models. Alternatively, doubly-charged Higgs 
(H ) may be produced in pairs in left-right symmetric 
models [9J. Under certain conditions, the H ±:iz primar- 
ily decays into two leptons, producing a final state of 
four leptons without $t or jets. In order to be sensitive 
to multiple new physics models, no $t or jet cuts are 
applied. The results of this search are interpreted us- 
ing an tip SUSY model within an mSUGRA framework 
only. A brief introduction to mSUGRA and lt p SUSY 
phenomenology is included below. 

In the minimal supersymmetric model [J], there are 
two Higgs doublets, as well as supersymmetric partners 
for every SM particle that differ by 1/2 unit of spin. The 
superpartners of the electroweak gauge bosons and Higgs 
bosons mix to produce four neutral mass states (x.1-4) 
and four charged states (Xi2)- The supergravity frame- 
work (mSUGRA) 0] has five free parameters: the univer- 
sal gaugino mass (M1/2), the universal scalar mass (Mo), 
the ratio of the Higgs vacuum expectation values (tan/3) , 
the trilinear coupling (Aq), and the sign of the higgsino 
mass parameter (/z). 

R-parity (R p ) is a quantum number defined such that 
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FIG. 1: Lowest order Feynman diagrams for Hp SUSY de- 
cay of the LSP via virtual sleptons governed by a) the A121 
coupling, and b) the A122 coupling. 



all SM particles have value +1 and all superpartners 
have value —1. Many searches assume that R p is con- 
served so that the LSP is stable, although there are vi- 
able SUSY models with Rp violation. One scenario in- 
cludes a sing le#p term in the renormalizable superpoten- 
tial, XijkLiLjEk, where Xijk are lepton number violating 
coupling strengths; and k are family indices ranging 
from 1-3; Li is the left-handed lepton doublet superfield; 
and Ei is the right-handed charged lepton singlet super- 
field. In this scenario, superpartners are produced in 
pairs at hadron colliders. Muon decay measurements Q 
constrain \i2k IS 0.14 so that heavier superpartners will 
cascade decay into the LSP with nearly 100% branch- 
ing ratio, producing at least two LSPs per event. The 
LSP, which is the Xi in mSUGRA for the values of M 1/2 
pertinent to this study, will then decay into two leptons 
plus a neutrino, as shown in Fig. [TJ resulting in a final 
state of four or more leptons. The detector acceptance 
limits our sensitivity to \\2k ft 3 • 10~ 3 , since a smaller 
value will result in the LSP decaying far from the in- 
teraction region At least two of the four leptons 
will be electrons (muons) for LSP decays governed by 
the A121 (A122) coupling. Final states with taus are not 
used since they are less efficiently identified. Recently, 
lower mass limits have been set by the D0 Collaboration 
on the Xi ranging from 115-119 GeV/c 2 and on the Xi 
from 229-234 GeV/c 2 in this type of itSUSY model, 
by searching for three leptons plus |£t uM- Previously, 
the limits from LEP were ~53 GeV/c 2 on the Xi an d 
103 GeV/c 2 on the X f 03 ■ 

The components of the CDF II detector relevant to this 
analysis are briefly described here; a more complete de- 
scription can be found elsewhere [l3| . The px and r\ [l4| 
of charged particles are measured by a silicon strip detec- 
tor [l5[ and a 96-layer drift chamber (COT) [l6| inside a 
1.4 T solenoidal magnetic field. The COT provides cover- 
age with high efficiency for |ry| < 1. For 1 < < 2, the 
silicon detector is mainly used. Electromagnetic (EM) 
and hadronic calorimeters surround the tracking system. 
They are segmented in a projective tower geometry and 
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measure energies of charged and neutral particles in the 
central (\r]\ < 1) and end-plug (1.1 < \rj\ < 3.6) regions. 
Each calorimeter has an EM shower profile detector po- 
sitioned at the shower maximum. Four layers of pla- 
nar drift chambers located outside the central hadron 
calorimeters and another set behind a 60 cm thick steel 
absorber detect muons with \rj\ < 0.6. Additional drift 
chambers and scintillation counters detect muons in the 
region 0.6 < \rj\ < 1.0. Gas Cherenkov counters [l7| 
measure the average number of inelastic pp collisions per 
bunch crossing and thereby determine the luminosity. 

Electron identification criteria include a reconstructed 
track associated with a cluster of energy in the EM 
calorimeter for I77I < 2.0. "Tight" central electrons and 
end-plug electrons are required to have high track quality, 
calorimeter cluster Et that is consistent with the track 
Pt, a high fraction (usually 95%) of the total calorimeter 
energy deposition in the EM section, and an EM lateral 
shower profile consistent with test beam measurements. 
A "loose" central electron sample is also collected by re- 
moving the shower profile requirements. Muons are iden- 
tified for 1 77 1 < 1 by a charged track consistent with being 
minimum ionizing in the calorimeters matched to a recon- 
structed track segment ( "stub" ) in one of the muon drift 
chambers. Stubless muons are also accepted for \r/\ < 1 if 
all other criteria are satisfied except for the stub require- 
ment and the track does not point to any of the muon 
drift chambers. All tracks associated with electrons or 
muons must be consistent with being produced at the 
pp collision point. Leptons are required to be separated 
from each other by AR = y/(Ai]) 2 + (A</>) 2 > 0.4 and 
to have \z\ < 60 cm at the point of closest approach to 
the primary vertex. Leptons are required to be isolated 
such that E^° divided by the E^ of the electron or p^ 
of the muon is less than 0.1, where E^° is the transverse 
energy within a cone of AR < 0.4 that is unassociated 
with the lepton. 

The data were collected with lepton triggers requir- 
ing at least one central electron (muon) candidate with 
E T > 18 GeV (p T > 18 GeV/c). We select events with 
two or more leptons that have £r (pt) > 20, 8, and 
5 GeV (GeV/c) for the first, second, and additional lep- 
tons, respectively. At least one identified lepton must 
pass the trigger requirements. From this sample, trilep- 
ton and >4-lepton events that pass additional cuts (de- 
scribed below) are used as signal candidates, while dilep- 
ton events are used to validate the background predic- 
tion. 

Several SM processes result in final states with three 
or more isolated leptons. The dominant contribu- 
tions include dileptons from the Drell-Yan Z/j* process 
(DY) with additional leptons from photon conversions 
or misidentified jets, and secondary contributions from 
diboson (WZ, ZZ) production. Except for those due 
to misidentified jets, the backgrounds are estimated us- 
ing Monte Carlo (MC) simulation. The simulated de- 
tector acceptances are not a 100% accurate representa- 
tion of the data and therefore are corrected by £data/eMC> 



where edata (cmc) is the efficiency for lepton reconstruc- 
tion, lepton identification, and photon conversion re- 
moval requirements measured in data (MC). This cor- 
rection varies from 0.6-1.0 per event, depending on the 
number and types of identified leptons. Lepton recon- 
struction and identification efficiencies are measured from 
DY and J/ip dilepton events in data and MC simulation. 
Trigger efficiencies are determined from W — ► ev events 
for electrons and from Z — > events for muons. The 
probability for a jet to be mis-identified as a lepton is 
determined as a function of Et for electrons and pt for 
muons from jet data samples, by counting the numbers of 
jets that remain after applying the lepton identification 
requirements. A similar procedure is described in more 
detail in [Hj]. 

Samples of fi p SUSY simulated events are gener- 
ated using pythia [nj within the mSUGRA framework 
with input of the supersymmetric couplings and particle 
masses from ISAJET [20j. For each point in the SUSY 
parameter space, the next-to-lcading-order cross section 
is calculated with PROSPIN02 [2l|. In all cases X2X1 
and xfxf production dominate and are the only pro- 
cesses generated. The Xi mass, which is roughly pro- 
portional to M1/2, has the largest impact on the signal 
acceptance. The x.i mass is also affected by the sign(/i), 
but has little dependence on Mo, tan/3, and A . There- 
fore, three of the five mSUGRA parameters are held con- 
stant: M = 250 GeV/c 2 , tan/3 = 5, and A = 0; while 
sign(/i) and Mi/ 2 are allowed to vary. Since the x.i mass 
decreases with larger values of Mq and tan/3, the chosen 
values are expected to produce conservative results. 

To reduce SM background contamination, an event is 
removed if any of the following conditions are met for 
opposite-sign, same-flavor lepton pairs: invariant mass 
between 76-106 GeV/c 2 ("Z veto"), invariant mass below 
15 GeV/c 2 , or 160° < A<j> < 200°. Identified cosmic ray 
events are removed. Electrons are removed if there is a 
partner track identified as coming from a photon conver- 
sion. An event is also rejected if the invariant mass of the 
two highest transverse energy leptons is below 20 GeV/c 2 
to reduce heavy flavor and radiative photon conversion 
backgrounds. In order to improve signal-to-background 
in the 3-lepton data sample only, we require that one of 
the two leading leptons be an electron when evaluating 
the A121 scenario, and a muon for the A122 scenario. This 
definition leads to a 22% overlap for the background and 
~15% overlap for the signal (depending on model pa- 
rameters) between the trilepton event samples. After all 
selection criteria, the signal acceptances are ~11% and 
~4% for the trilepton and >4-lepton data samples, re- 
spectively. 

Uncertainties are determined separately for SM back- 
ground and $p SUSY expectations, in each data sam- 
ple. The sources of systematic uncertainty on the trilep- 
ton background prediction include jets misidentified as 
leptons (12-13%), lepton identification efficiency (5-6%), 
luminosity measurement (6%), choice of parton distribu- 
tion functions (2%), cross section (5-6%), photon con- 



6 




e T (GeV) 

FIG. 2: ]Et for events in the trilepton signal samples. For 
comparison, an ft p SUSY sample is shown with M(xi) = 99 
GeV/c 2 . 



version identification (3%), and initial state radiation 
(4%). For events with >4-leptons, the uncertainty on the 
background is dominated by jets misidentified as leptons 
(41%) and Z/ 7 * + 7 MC statistics (38%). The total sys- 
tematic uncertainty on the signal acceptance ranges from 
12-15%, based on the number and types of leptons, where 
the largest contribution is due to the uncertainty on the 
low-pT (< 20 GeV/c) lepton identification efficiency. 

The background prediction is validated through the 
use of control regions in which one or more of the event 
selection criteria are inverted. For each control region, 
the SM prediction is compared to the data, as shown in 
Table [T] In the ee and dilepton control regions, the 
background is dominated by DY. In the case of e\x con- 
trol regions, the largest background is Z/j* — > tt; how- 
ever, there are also significant contributions from WW, 
it, and jets misidentified as leptons. For trilepton con- 
trol regions, the largest backgrounds originate from DY 
events where the third lepton is due to either a photon 
conversion or a misidentified jet. The agreement between 
predicted and observed events in the control regions indi- 
cates that the backgrounds are validated, and we proceed 
to examine the signal data samples. 

In the trilepton data samples, a total of 6 events 
are observed: 5 events with an expected background 
of 3.1 ± 0.7 (stat.) ± 0.4 (syst.) events for the Am 
scenario, and one event with an expected background 
of 1.9 ± 1.0 (stat.) ± 0.3 (syst.) events for the 
A122 scenario. The distribution of the observed 
events are consistent with the background prediction, 
as shown in Fig. [2j No events are observed in the 
>4-lepton data sample, with an expected background of 
0.008 ± 0.003 (stat.) ± 0.003 (syst.) events. We interpret 
the results as being consistent with the hypothesis of no 
signal and therefore set mass limits on the x? aud Xi par- 
ticles. The cross section limits are calculated by com- 
bining the 3-lepton and >4-lepton signal samples using a 



TABLE I: Predicted and observed events found in the control 
and signal samples. The leptons are listed in order of decreas- 
ing transverse energy. The first lepton must either be a tight 
central electron or "stub muon." Systematic and statistical 
uncertainties are added in quadrature. 



Region Criteria Failed Predicted Observed 

Background Events 

Two-Lepton Control Samples 

ee Z veto 13948 ± 1536 14019 

ee 2142 ± 230 2125 

llll Z veto 7474 ± 809 7499 

llll 1264 ± 141 1339 

e/i 117.6 ±12.9 112 

ae 186.8 ± 22.9 203 



Three-Lepton Control Samples 



eel 


Z veto 


8.8 ± 1.9 


12 


Llfll 


Z veto 


5.1 ± 1.2 


2 


efil 


Z veto 


0.55 ±0.04 





III 


Z veto 


14.4 ±2.9 


14 


eel 


A</> only 


2.1 ± 0.3 


2 


liliI 


Acf> only 


1.2 ±0.2 


4 


efil 


A(f> only 


0.35 ±0.04 





III 


A<f> only 


3.7 ±0.3 


6 



Four-or-More-Lepton Control Samples 

llll Z veto 0.15 ±0.02 

llll A(j> only 0.006 ±0.003 

Three-Lepton Signal Samples 

A121 scenario 3.1 ± 0.8 5 

A122 scenario 1.9 ± 1.0 1 

Four-or-More-Lepton Signal Sample 

A121 , A122 scenarios 0.008 ± 0.004 



TABLE II: Expected and observed 95% C.L. limits on the Xi 
and xt masses for 3- 10~ 3 < Ai 2fc < 0.14, Mo = 250 GeV/c 2 , 
tan/3 = 5, and Aq = 0. 



ft p SUSY 


M(x?) [GeV/c 2 ] 


M(xr 


) [GeV/e 2 ] 


scenario 


exp. 


obs. 


exp. 


obs. 


Am, n > 


105 


102 


192 


185 


A121, £i < 


101 


98 


192 


186 


A122, M > 


108 


110 


198 


203 


A122, £t < 


103 


106 


195 


202 



multi-channel Bayesian method similar to [22| , shown for 
one A p SUSY scenario in Fig. [3l The limit calculation 
takes into account correlations between the 3-lepton and 
>4-lepton samples. The resulting mass limits, at 95% 
confidence level (C.L.), are presented in Table HT1 

In conclusion, we have performed a search for anoma- 
lous production of events with three or more leptons us- 
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FIG. 3: Experimental 95% C.L. cross section limits and theo- 
retical ftp SUSY cross sections versus LSP mass, for 3- 10 -3 < 
A122 < 0.14, M = 250 GeV/c 2 , tan/3 = 5, A = 0, and /Lt > 0. 

ing a sample of CDF Run II data corresponding to 346 
pb _1 of integrated luminosity. Finding results consistent 
with the SM, we set limits on the x? and X\ masses 
within an mSUGRA framework. The \i mass limits 
range from 98 to 110 GeV/c 2 , while the chargino mass 
limits range from 185 to 203 GeV/c 2 at 95% C.L. de- 



pending on the choice of model parameters. These re- 
sults significantly improve upon the LEP limits [l2j and 
are comparable to the D0 limits [ll|. While no evidence 
for new physics was found, the >4-lepton data sample is 
virtually background free and provides an excellent tech- 
nique for detecting new physics with more luminosity in 
the future. 
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